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A simplified test of universality in Lattice QCD is performed by analytically evaluating the contin- 
uous Euclidean time limits of various lattice fermion determinants, both with and without a Wilson 
term to lift the fermion doubling on the Euclidean time axis, and comparing them with each other 
and with zeta-regularised fermion determinant in the continuous time — lattice space setting. The 
determinant relations expected from universality considerations are found to be violated by a certain 
gauge field-dependent factor, i.e. we uncover a "universality anomaly". The physical significance, 
or lack thereof, of this factor is a delicate question which remains to be settled. 
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The low-lying spectrum of the naive lattice Dirac operator approximates the low-lying spectrum of the continuum 
Dirac operator but with a 16-fold degeneracy due to fermion doubling [1]. Lattice QCD with a naive fermion is 
r T | ' therefore regarded as a regularisation of continuum QCD with 16 degenerate fermion flavours. Lattice QCD with a 
staggered fermion [2] is regarded as a regularisation of 4 flavour QCD, in accordance with the fact that one naive 
fermion flavour is equivalent after spin diagonalisation to 4 staggered fermion flavours [3] . Lattice QCD with a Wilson 
fermion [1], where the Wilson term is added to the naive fermion action to lift the fermion doubling, is regarded as 
a regularisation of QCD with a single fermion flavour. Implicit in this is a universality hypothesis: the LQCD's with 
naive fermion, staggered fermion, and Wilson fermion are all in the right universality class to reproduce continuum 
QCD, the only difference being in the number of continuum fermion flavours the different lattice fermion formulations 
describe. 
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' It is highly desirable to test this universality hypothesis wherever possible. This is particularly important in view 
of the fact that LQCD calculations with both dynamical Wilson and staggered fermions are currently being pursued 
at great effort and expense (see [4] for a recent review with references to the literature). An interesting quantity to 
consider in this context is the fermion determinant, which appears in the QCD functional integral when the fermions 
are dynamical. In LQCD with dynamical staggered fermions, the fourth root of the staggered fermion determinant 
is used as the fermion determinant for a single quark flavour. An important test of the universality hypothesis would 
therefore be to check whether the fourth power of the Wilson fermion determinant coincides with the staggered fermion 
determinant in the continuum limit, or, equivalently, whether the 16th power of the Wilson determinant coincides 
' with the naive fermion determinant in this limit. Such a test appears analytically impossible with currently known 
, techniques though. However, a simplified version of this test is feasible: instead of the full continuum limit one can take 
the continuum limit for one of the spacetime coordinates while keeping the remaining coordinates discrete. The effect 
of including in the fermion action a Wilson-type term to lift the fermion doubling along the chosen coordinate axis 
can then be investigated. In this paper we perform such a simplified test of universality by evaluating the continuous 
Euclidean time limits of the various lattice fermion determinants, both with and without the time part of the Wilson 
term in the action, and comparing them with each other and with the fermion determinant in the continuous time — 
lattice space setting. (The fermion determinant in the latter setting is defined via zeta-regularisation.) 

On a finite spacetime lattice, with N a sites along the spacial axes a = 1, 2, 3 , Np sites along the Euclidean time 
axis, time lattice spacing=a, spacial lattice spacing=a', setting (3 = aNp— time length (=inverse temperature in the 
finite temperature QCD setting), we consider lattice fermion actions of the form S = a(a') 3 r % $(x, r)D^ r '^(x, r) 
where (x, r) runs over the lattice sites and 

D[r) = 74 ±V 4 + ^A 4 + D space + m (1) 

V 4 ^(x, t) = |(tT 4 (x, t)?/>(x, t + a) - U 4 (x, t - a)- 1 ^>(x, r - a)) (2) 

A4^)(x, t) = 2ip(x, t) — U±{x, t)i/j(x, t + a) — U^x, r — a)~ ip(x, t — a) (3) 

n( r ') - ±yf + sL A (A) 

J - / space ~ a'/* space ~ 2a ' ^space 
3 

rfspacetp(x,T) = ^ 7 CT \ (U a (x, t)V>(x + <7, t) - U a {x - <T, t) ~ 1 lj}{x - <T,t)) (5) 
<y=l 
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A space -0(x,T) = ^(2-0(x,t) - J7 ct (x,t)t/;(x + (t,t) - J7 ct (x-<t,t) 1 ^-g,t)) (6) 

cr=l 

For r = r' = this is the naive fermion action while for r = r' ^ it is the Wilson action. We are going to evaluate the 
continuous time limit (a — > , Np — » oo with /3 = aATg held fixed) of the fermion determinants detD^ 1 at r = and 
r = 1 and compare them with each other and with the zeta- regularised fermion determinant det^D a in the continuous 
time — lattice space setting. The Dirac operator in the latter setting is 

D = 74 (£ + Ai) + D space + m (7) 

with A4(x,t) being the 4-component of a smooth continuum gauge field such that {/ 4 (x, t) = Te^o °- A 4( x > T +( 1 *) a ) dt 
(T=t-ordering) is the lattice transcript. The subscripts "a" in £>i r ^ and D a refer to the operators defined by replacing 
U4 —> e~ aa U~4 in (l)-(6) and A4 — > A 4 — a in (7), respectively. The role of the complex parameter a is to incorporate 
the effect of a general boundary condition at the time boundaries: Da^ (resp. D a ) with periodic time b.c. has the same 
spectrum and determinant as D^ r ' (resp. D) with time b.c. ^>(x, /3) = e a ^ip(x.,0). Thus the introduction of a allows 
us to always take periodic time b.c. when considering the fermion determinant. It can also be used to incorporate a 
chemical potential fi : QCD at finite temperature and density, where the fermion fields satisfy anti-periodic time b.c, 
corresponds to a = fi + in/ (3. The gauge fields are required to satisfy periodic time b.c. The spacial b.c.'s for the 
fermion and gauge fields will not play a role in our considerations; we can simply take them to be periodic. 

The term ^A 4 in (1) is the "time part" of the usual Wilson term. It lifts the fermion doubling on the Euclidean 
time axis when r =/= 0. Therefore, if we think of the continuous time — lattice space setting as the "continuum setting" , 
then the aforementioned universality hypothesis, which relates the continuum limits of the naive, staggered and Wilson 
fermion determinants, translates into the following 
Simplified universality hypothesis: 

lim detD^ = ( lim detD^A (mod p.i.f.'s) (8) 

where "p.i.f.'s" refers to 'physically inconsequential factors'. This is now something which can be checked analytically. 
Our main technical result in this paper is the following 
Result of calculation: 

lim detD^ = ( lim deiD™)* e~ fa Tr & A -*"' (T))dT (mod p.i.f.'s) (9) 

where A space (r) is defined on the space of lattice spinor fields V( x )> living only on the spacial lattice, by replacing 
^(x, t) by V'(x) m (6). The p.i.f.'s in (9) are gauge field-independent factors whose only effects are to produce constant 
(vacuum) shifts in certain physical quantities. They include inverse powers of a which diverge in the a — > limit. 

The result (9) reveals a "universality anomaly" : the exponential factor in the right-hand side violates the universality 
expectation (8). Thus it is important to ascertain the significance, or lack thereof, of this factor. Since it is gauge field- 
dependent it cannot strictly speaking be regarded as a p.i.f. in the continuous time — lattice space theory. However, 
since the spacial Wilson term ^p-A space formally vanishes in the spacial continuum limit, one could argue that the 
exponential factor is effectively a p.i.f. when one goes on to take that limit. This is a delicate issue though, since 
Tr -^-Aspace actually diverges in in the a' — > limit (note that in the free field case the largest eigenvalue of jrA space 
is ~ i). Further study is required to clarify this issue. 

Since the reasoning which leads to the universality expectation (8) is the same as that which leads to the expectation 
that LQCD with naive, staggered and Wilson fermions are all in the same universality class, (9) is a potential reason for 
concern about whether the latter universality hypothesis really holds. It would therefore be a significant reassurance 
if the anomaly factor in (9) can be shown to be physically inconsequential. It should be noted, however, that even 
if this turns out not to be the case, it would not in itself invalidate the universality hypothesis for LQCD with 
naive, staggered and Wilson fermions, since the comparison between naive fermion (r = r' = 0) and Wilson fermion 
(r = r' ^ 0) is not covered by (9) and we could still be lucky in this case. But it would certainly raise a serious 
concern. 

Remarkably, the anomaly in (9) is mirrored by an ambiguity in the zeta-regularised fermion determinant in the 
continuous time — lattice space setting. The latter can be expressed either as det^D a or det^^Da) (since r/S = ^*74)- 
Formally, the determinants of D a and 7 4 -D Q coincide, but it turns out that the rigorously defined zeta-determinants 

do not. In fact we find det ( (j 4 D a ) = e - ' f» Tr ( ^ rAspa ^ (T)) dT det ( D a (mod p.i.f.'s), and 
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lim detD^ = det c D a , lim detD^ = detr{^D a f (mod p.i.f.'s) (10) 

This shows that the lattice regularisations are consistent with zeta-regularisations of the fermion determinant in 
the continuous time — lattice space setting, and that the requirement that the anomaly factor in (9) be physically 
inconsequential is also necessary for consistency of continuous time-lattice space QCD when the fermion determinant 
is defined by zeta-regularisation. 

In the remainder of this paper we sketch the derivation of (9)-(10) and give other, more explicit, expressions for 

the a — > limits of detDa^ and detDa\ The full details are provided in [5]. It is convenient to regard ip(~x., r) as a 
function ^(r) living on the lattice sites of the Euclidean time axis and taking values in the vectorspace W = {^(x)}, 
i.e. the space of lattice spinor fields living on the spacial lattice only. Set N := dimVF. Define the linear operator 
U±{t) on W by ([/4(t)V>)( x ) = U^x, r)^(x). The operator D space (T) on W is defined similarly by replacing ip(x, r) by 
V>(x) in (4)-(6). Since *(/3) = *(0) we can represent <f by the vector * = (*(0), . . . , $(Np - 1)) where *(fc) = ^(ka). 
Then is represented by 

flM*(ife) = d£l(k)$(k - 1) + d { r \k)y(k) + d[ r \k)^(k + 1) (11) 

where the operators d^ r \k) :W^W are given by d^\k) — ^- (74 - r) Ui(k) , d^\(k) = (74 + r) J7 4 (fc - 1) _1 , 
d!^\k) = - + D S pace{k) + to with Ui(k) := Ui(ka) and D S pace(k) := D space (ka). The generalisation of Z)( r ) to Z)i r '', 

given by U4 — > e~ Qa [/4 , is equivalent to rf^ — ► e^" ^} in (11). After writing as an Np x iVjg matrix, its 
determinant can be straightforwardly evaluated via the method of [6]. The cases r = 1 and r ^ 1 require separate 
treatments due to the fact that d±\(k) is invertible when r ^ 1 but not when r = 1. The details of the calculation 
are provided in [5]; here we simply quote the results, assuming for convenience that Np is even in the 1 case: 

detD™ = {^) NN " det{[\ J)-^*^)) (12) 



det 



= J e -«/^/2 J-j" det {x + aM{k)) det(l-e af) V(Np)) (13) 

where M(k) := ^A(fc) + to (i.e. the scalar part of D S pace(k) + to). The linear maps U^{Np/2) on W(BW and 
V(Np) on 14 7 in (12)-(13) are defined as follows. Exploiting the periodicity of the link variables to define dj(k) for all 
k S Z, periodic under fc — > fc + iVjg , and thereby define D^ r ^(k) for all fc £ Z , we consider the equation D^^!(k) = 
(no periodicity requirement on $(k)). In the r / 1 case, since the <4i( fc )' s 

are invertible, it is clear from (11) 

that solutions ^>(k) are specified by two initial values. Thus the solution space is isomorphic to W © W. Setting 
$i(n) = i&(2n) and ^>2{n) = *fr(2n+l) the solutions are determined from their initial values via an evolution operator: 

(^(n)) = ^ r \ n ) (^(o)) " ^ e °P era t° r (Np/2) appearing in (12) can also be characterised as follows: Due to 

the Ayi-periodicity of the d^(fc)'s in (11) there is a linear map on the solution space defined by ^(k) 1— ► ^>(k + Np), 

or, equivalent^, (*i(n), * 2 (n)) >-> (*i(n + Np/2),$ 2 {n + Np/2)). This map coincides with U^{Np/2) when the 
solution space is identified with W © VF. 

In the r = 1 case, solutions to D^^(k) — are determined by just a single initial value; this is connected with 
the non-invertibility of d±\(k) and can be seen, e.g., from the expression (16) below. Thus the solution space in 
this case is isomorphic to W. The evolution operator V(k) determines solutions from their initial value through 
if?(k) = V(k)^>(0). The V(Np) in (13) can be alternatively characterised as the linear map on the solution space which 
maps if(k) i-> 4?(k + Np). 

Finite difference approximations to differential operators in one variable and their determinants have been studied 
in [7,8] and we are going to use a convergence result from there. In the setting of [7,8], specialising to 1st order 
differential operator, the operator L and its finite difference approximation L have the forms 

L = L^t)-^ + L (t) , L = L 1 (k)-d + L (k) (14) 
or a 

(r e R, k e Z) d e {d+,d-} , 0+tf(fc) = *(fc + l) - tf(Jfe), 0-*(fc) = tf(fc) - tf(fc-l), with L^L^k) :W^W 
being periodic under r^r + /3, k — > k + Np , respectively, and 
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Lj(k) = Lj(ka) + 0(a) (j = 0, 1). 



(15) 



Then the solutions to L^f(r) = and L^>(k) — are both determined by a single initial value, so the solution spaces 
in both cases are isomorphic to W. Solutions ^> approximate solutions $ , i.e. if "&(0) = ^(0) then ^>(k) « ^>(ka) for 
small a. Consequently the evolution operator U(k) for L^f — approximates the evolution operator U(j) for L'fy = 0. 
(Explicitly, U(t) = Te~f» L 1 (t)- 1 L (t)dt^ ^ particular onc has thc following (cf. §3 of [8]): 
Convergence Theorem: U(Np) — > U((3) for a — > with aiV/3 = /3 held fixed. 

An obvious variant of this which we will make use of is the following. If p is a multiple of Np and L = L\ ^ d + L 

with the Lj(k)'s periodic under k — > k + Np/p and satisfying Lj(k) = Lj(kpa) + 0(a) (j = 0, 1) then U(Np/p) — > W(/3) 
for a — > 0. Furthermore, if W is replaced by W\ © W2 in the preceding then the convergence theorem continues to 
hold when d is replaced by with d, d G {d + , d~}. 

In order to apply the convergence theorem to evaluate the a — > limits of (12)-(13) we need to rewrite Z)^ r - ) in the 
form of L in (14), or its aforementioned variant. We have only been able to do this in the r = and r = 1 cases. The 

(r) 

problem of evaluating lim a ^o detD a in thc general r case therefore remains for future work; new techniques beyond 
those of [7,8] may be required for this. In thc r = 1 case we specialise to a 7-representation where 74 = ^ J J 3 ^ and 
decompose TU = W+ TU- so that 74 = ±1 on W±. Then, in terms of this decomposition, 



a 



tf>«-*r' (( *; iw ~\ ° 1 (it) 



4 1, W = 74( " (1 - Cf4(( *- 1) ° r ' ) , ° JtD^Mt-. (18) 



where 

V U 4 {ka) / 

(k-l)a)- 1 ) 

I (f / 4(fca )-l) 
Clearly L*- (fc) is periodic under k — > fc + iVyg and 

Lp ) (fe)= 74 + O(a) , 4 X) (fc) =7 4 A 4 (fca) + + ™ + O(a) ■ ( 19 ) 

The convergence theorem now gives lim a ^ V(Np) = V(/3) where V(r), acting on W, is the evolution operator for 
D^>(t) = with D being thc Dirac operator (7) of the continuous time — lattice space setting. Using this and noting 
det(l + aM(k)) = e aTrAl ^ + 0(a 2 ) thc a -> limit of (13) is now obtained: 



lim deW^ = ( l -) NN * e -«W2+if:TrM { r )d r ^ ^ (2Q) 



• 1 \ NN 

The gauge field-independent factor l/a NNf3 , which diverges in the a — * limit, is physically inconsequential; it can at 
most give rise to an overall constant shift in the calculation of certain physical quantities (such as the energy density 
in finite temperature QCD). 

An application of the zeta-regularised determinant formula for differential operators in one variable, Theorem 1 of 
[7], leads to an expression for det^D a which coincides with (20) without the l/a NNf> factor and with M(r) replaced 
by ±M(t) (the details of this are given in [5]). The sign ± depends on the choice of cut in the complex plane used 
to define thc zeta-determinant. Choosing this so that the sign is "+" wc then have lim Q ^ a NN/3 detD^ = det^D a 
which establishes the first part of (10). 

In thc r = case, with ^ represented by (^i(n), ^(rc)) as before, we have 



1 2a I d- ' 



(21) 



where 
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m (n \-( ° ]l{n) \ f{0) (D space (2na) + m K ( n ) \ 

Ji(n) =74t7 4 ((2n- l)a) _1 , ifi(n) = 74 J7 4 ((2n + l)a) (23) 
JoW=74^(t/4(2na)-[/ 4 ((2n-l)a)- 1 ) , ^ (n) - 74 ^ ( ^ 4 ((2n + l)o) - t/^na)" 1 ) (24) 

Clearly L^\n) is periodic under n — > n + Np/2 and 

1 V 7 V74 Oy W V y V 74^ 4 (2na) D spoce (2na) + m J W V 7 

The convergence theorem then gives lim a _> U(Np/2) = U([3) where U{(3), acting on W® W, is the evolution operator 
for 5(£)(r)=0 with 

5= ( D space { T )+m 74 (A + A 4 (t))\ 
V74(^ + A 4 (t)) fl spm W+™/' 

Introducing O = (_\ on W © W, we find after a little calculation O' 1 D O = ( (t*^)" 1 ^ (7475) « ^ w h e re D is 
the Dirac operator (7). It follows that U(J3) = C" 1 ( (t 4 ^)" 1 (7475) a Using thiS; the a ^ limit of (12) 

is now obtained: 

lim detD^ = ( 1_) NN? e -<*PN det 2^ _ e «/3 V (/3)) (27) 
a^o V 2a/ 

Again there is a physically inconsequential, divergent factor, (l/2a) NN P . Comparing (27) with (20), and noting that 
e'l TV m(t) dr = emf3N/2 Jj a Tr{^r\ spacc (r))dr where is & p j ^ we obtain thc claimed result (g) Further- 

more, an application of Theorem 1 of [7] gives det^j^Da) = e"! 1 ? 1 )! 3 " 1 / 2 det(l — e Q/3 V(/3)) where the sign =p again 
depends on the choice of cut in the complex plane [5]. It follows that lim a ^o {2a) NNf3 detD^ = e ±l3aN dei<j(7 4 D Q ) 2 . 
Since e ±/3aJV is a p.i.f. this establishes the second part of (10). 

In thc free fcrmion case, Z — detDa^ with a = /1 + iir//3 is the partition function of a free lattice Dirac fermion 
gas at temperature T = 1/(3 and chemical potential fi. After noting det(l — e al3 V(/3)) = det(l + e~^ H ) where 
H = 7 4 (D space + m) it is easy to check that physical quantities such as the mean energy (E) = --§plogZ (with /3jt* 
held constant) calculated from the expressions in this paper are consistent with the standard results described, e.g., 
in [9]. One subtlety is that in thc r = 1 case lim a ^o (E) has an additional shift of — 2 ^ p (^p-A space (p) + m) due 
to the exponential factor in (20), while such a shift does not occur in thc r = case (cf. (27)). This shift can also 
be seen from the standard calculation in §18.11 of [9]: There the computation is reduced to a weighted sum of the 
residues of a certain function g(z)/z. In addition to the poles at z — z± discussed there, there is also a pole at z = 
and its residue is readily seen to give the aforementioned shift in lim a ^ (E). In the r = case there is no pole at 
z = and hence no corrsponding shift in that case. 

We have considered only the simplest lattice fermion formulations in this paper, namely the naive and Wilson 
fermions (with the connection to staggered fermions due to thc relation between the naive and staggered fermion 
determinants). More recent formulations of major current interest are the domain wall [10] and overlap [11] fermions, 
as well as other chirally improved approaches to lattice fermions [13]. It appears feasible to evaluate the continuous 
time limit of the domain wall fermion determinant (with finite, discretised fifth dimension) along the same lines as 
here; this is currently under investigation. Doing this in the case of overlap fermions seems more technically difficult 
though. We remark that the overlap fermion determinant has previously been successfully tested (in both the chiral 
and vector cases) against known continuum expressions in 2 dimensions [12]. 
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